Half Metallicity in Hybrid BCN Nanoribbons 
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We report a first-principles electronic-structure calculation on C and BN hybrid zigzag nanorib- 
bons. We find that half- metallicity can arise in the hybrid nanoribbons even though stand-alone C 
or BN nanoribbon possesses a finite band gap. This unexpected half- metallicity in the hybrid nanos- 
tructures stems from a competition between the charge and spin polarizations, as well as from the 
7r orbital hybridization between C and BN. Our results point out a possibility of making spintronic 
devices solely based on nanoribbons and a new way of designing metal- free half metals. 
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The discovery of low-dimensional carbon allotropes 
such as fullerenes [l( and nanotubes [2[ has stimulated in- 
tensive research on metal-free magnetism, owing to their 
small spin-orbit coupling and long spin scattering length. 
Since the first report of room-temperature weak ferro- 
magnetism in polymerized Ceo 0, magnetism in pure 
carbon materials has been demonstrated by several ex- 
perimental groups [1, H, @] . Carbon defects or adatoms 
are believed to be the origin for the unexpected mag- 
netism 0,H@|. 

A question arises that can low-dimensional carbon 
nanostructures exhibit half metallicity? Half metals are 
ideally for spintronic applications because they have one 
metallic spin channel and one semiconducting or insu- 
lating spin channel. Techniques for tuning electronic 
properties of graphene nanoribbon (GNR) have been ad- 
vanced [To|, EH EH > for example, by terminating a single 
graphite layer in one direction. Recently, Son et al [13[ 
predicted, based on density functional theory (DFT) cal- 
culation, that zigzag edge GNRs (ZGNRs) can be con- 
verted to half metal when an external transverse elec- 
tric field is applied onto the graphene layer. Our recent 
study using hybrid DFT confirmed their prediction [14j. 
However, we found that a very strong field is required to 
achieve half metallicity, a limitation for wide application. 

In the absence of external field, ZGNRs are semicon- 
ductor with two localized electronic edge states [Til, EH, 
17l flil . [Tit - These two ferromagnetically ordered edge 



states are antiferromagnetically coupled. In the presence 
of an external field, an electrostatic potential difference 
is generated between the two ZGNR edges, which causes 
band crossing and charge transfer between the two edge 
states. As a result, the ZGNRs are converted to half 
metal from semiconductor. Since the critical step to 
make half-metallic ZGNR is to break the symmetry of 
the two edge states, an alternative approach is to use 
chemical decoration. For nanoelectronic applications, the 
chemical approach has many advantages over the strong- 
field approach. 



In this Letter, we present two nanostructure designs to 
achieve half-metallicity for carbon nanoribbons without 
using metal dopants. The first design is by implanting 
a BN row into a ZGNR. The principle of this materials 
design is based on the following observation: both the 
highest valance band (VB) and the lowest conduction 
band (CB) of ZGNR are mainly contributed by the two 
edge states which have opposite spin but are degenerate 
in energy. In contrast, for the zigzag BN nanoribbon, the 
highest VB is originated from the 7r orbitals of the N edge 
atoms, while the lowest CB is mainly localized at the B 
edge atoms. These results suggest that the symmetry of 
the ZGNR band structure near the Fermi energy can be 
broken by introducing a segment of BN nanoribbon. 

A geometry of the newly designed hybrid C/BN 
nanoribbon is shown in Fig[T^, where the armchair atom 
rows are normal to the ribbon direction while the zigzag 
atom chains are in the ribbon direction. At the two edges, 
H atoms are added to passivate the dangling a bonds. 
The hybrid ribbon structures are denoted as n-C^BN, 
where i refers to the number of armchair C rows in a 
unit cell, and n is the width of the nanoribbon in term 
of the number of zigzag chains. On one edge, B is the 
outmost non-hydrogen atom; this edge is denoted as the 
B-edge. The adjacent C atoms to the B-edge are named 
Cb- Similarly, the C atoms adjacent to the N-edge are 
named Cat. 

The DFT calculations were carried out using the VASP 
package [2pj ]. The projector augmented wave (PAW) 
[2l| method in Kresse-Joubert implementation [22] was 
used to describe the electron-ion interaction. The plane 
wave cut-off energy was set to 400.0 eV. The convergence 
thresholds for energy and force were 10 -5 eV and 0.01 
eV/A, respectively. Perdew-Wang functional [23[ known 
as PW91 was used in the generalized gradient approxima- 
tion (GGA). In a previous study of ZGNRs, we showed 
that the GGA gives qualitatively the same results as hy- 
brid density functional [14|. Here, the calculated band 
gaps of the pristine ZGNR and BN nanoribbon with 
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8 zigzag atom chains are 0.43 and 4.2 eV respectively, 
in_good agreement with previous theoretical calculations 

The calculated electronic structures of 8-C^BN (i=l, 
2, 3) are shown in Fig 1. The band structures of 8-CiBN 
and 8-C3BN exhibit different patterns compared to 8- 
C2BN, which can be understood from the Brillouin-zone 
folding. When i is even, the largest gap locates at T 
point, and when i is odd the smallest gap locates at T 
point. 8-C1BN has a nonmagnetic ground state with a 
small band gap of -0.1 eV. Both 8-C 2 BN and 8-C3BN 
are spin-polarized, and the latter is half metal since its 
valence and conductive bands overlap at the Fermi level 
in the spin-down channel. This half-metallicity is unex- 
pected, considering that both C and BN zigzag nanorib- 
bons possess a finite band gap. We further examined 
nanoribbons with different widths (n = 6 — 14). As shown 
in Fig Id, all QBN (i=l, 2, 3) nanoribbons become half 
metal when n > 12. The critical widths are n =12, 10, 
and 8 for z=l, 2, and 3, respectively. 
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FIG. 1: Band structures for (a) CiBN, (b) C 2 BN, (c) C 3 BN 
with 8 zigzag chains. The red (blue) lines represent the spin- 
up (spin down) channel, (d) The spin-resolved band gaps of 
n-CiBN nanoribbons. The filled square and hollow circle de- 
notes the spin-up and spin-down channel, respectively. (e)The 
atomic configuration of C^BN ribbons. Green, pink, grey, and 
sapphire balls denote carbon, boron, nitrogen, and hydrogen 
atoms, respectively. 

The calculated partial density of states (PDOS) indi- 
cates that the 8-C1BN has no spin polarization. Near 
the Fermi level, the occupied states are contributed by 
Cat and N edge atoms, while the unoccupied states are 
contributed by and B edge atoms. This is a charge 



polarized state with electron transfer from Cb to Cat. 
The charge polarization is clearly seen in Fig. 2a-g; the 
partial charge density from 0.05 eV below the Fermi level 
to the Fermi level is localized at the N-edge. 

In pristine ZGNR, every edge C atom at the edge has 
only two C neighbors and thus 1/3 nonbonded extra tt 
electron, which gives rise to the edge states and a high 
electron density at the Fermi level. A stablization mecha- 
nism is required to have a lower-energy ground state [13] • 
In the charge-polarization mechanism, the extra electron 
transfers completely from one edge to the other. In the 
spin-polarization mechanism, the C atoms at the two 
opposing edges have spin-opposite electrons, which low- 
ers the on-site Coulomb interaction U between electrons 
with opposite spin at the same site. The spin mechanism 
is more favorable in energy [17[. However, the charge- 
polarized state can be further stabilized by an external 
electric field. Hence, a ZGNR can be changed from hav- 
ing spin polarized state to charge polarized state under 
external field. In either the spin or charge polarized state, 
a ZGNR is semiconductor. Here, the half-metallicity for 
ZGNR is found in the region where the two states com- 
pete with each other 

The charge-ordered ground state of 8-C1BN is sta- 
blized by the mixing of tt orbitals of C and BN. In CiBN, 
each Cb atom is covalently bonded to a N atom, while 
each Cat atom is connected to a B atom. As shown in 
Fig. [2^, the occupied tt orbital of N is much lower in en- 
ergy than the unoccupied tt orbital of B. The tt orbitals 
of C locate in the middle. The interaction between the 
tt orbitals of C^ and N results in a higher antibonding 
orbital. However, the electron tends to occupy the lower 
bonding orbital between Cat and B. Such an orbital in- 
teraction leads to electron transfer from C^ to Cat. 

As the nanoribbon width (n) increases, the Coulomb 
interaction between the two edge states decreases, and 
thus the stabilization energy of the charge-polarized state 
also decreases. On the other hand, the on-site Coulomb 
interaction is a local property, and the spin-polarization 



energy is independent of the ribbon width [17[. Hence, 
the spin polarization becomes relatively more competa- 
tive as the nanoribbon width increases. 

As shown in Fig [2)3, when the ribbon width is not wide 
enough (e.g., n = 10), the ribbon still has an occupied 
N-edge and nearly empty B-edge state. However, the 
spin degeneracy is already broken. The spin-up and spin- 
down electrons tend to occupy one of the two subsites of 
the honeycomb lattice, which leads to the energy of spin- 
down edge state higher than of the spin- up one. When 
n=12, the occupied spin-down state at the N-edge and 
the unoccupied spin-down state at the B-edge cross over 
in energy. The B-edge is partly occupied, which weak- 
ens the charge polarization. More importantly, in this 
case, the spin-down channel becomes metallic. There- 
fore, the competition between the spin and charge po- 
larization leads to half-metallicty in this hybrid-ribbon 
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structure. As the ribbon width increases, a decrease of 
charge polarization (Fig. ^-c) and an increase of spin 
polarization (Fig. [2]i-f) can be clearly seen. 

For QBN with i > 1, the ir orbital hybridization only 
occurs between BN and its nearest neighboring C atom 
row. In the i=oo limit, the electronic structure of C^BN 
will converge to that of ZGNR. For different z, we ana- 
lyzed the competition between charge and spin polariza- 
tion and the half-metallicity based on the spin-polarized 
ZGNR ground state (see Fig. Eh). In ZGNR, C atoms 
at the two edges are degenerate in energy but with oppo- 
site spin. When the BN row is introduced, the energy of 
the bonding tt orbital between Cat and B is lower than 
that of the edge state in ZGNR. This bonding state is 
occupied by spin-up electrons. On the other hand, for 
Cb and N, both the bonding and anti-bonding states are 
occupied by the spin-down electrons. Hence, in QBN 
without charge polarization, the valence band maximum 
(VBM) stems from the anti-bonding states between N 
and Cb in the spin-down channel, while the conduct band 
minimum (CBM) is composed of the bonding states be- 
tween B and Cat in the same spin channel. With the 
increase of the ribbon width, the ZGNR band gap is re- 
duced, and the VBM and CBM are closer in location. 
Once the VB and CB overlaps at a specific point (i=3 
and n=8), spin-down electron at the B-edge will transfer 
to N-edge. Consequently, at both edges, the spin-down 
channel is partially occupied, and the partial charge po- 
larized nanoribbon becomes half metal. 

Note that implantation of a BN row in carbon nan- 
otube has been realized in the laboratory [H, [2fj, but 
not yet in ZGNR. Our first design for the hybrid C/BN 
nanoribbon is a highly ordered nanostructure. It may be 
challenging to realize this nanostructure in the labora- 




FIG. 2: The partial charge density within energy range [E/- 
0.05, E f ] in eV for n-CiBN with n equals to (a) 8, (b) 10, and 
(c) 12. Isovalue is 0.02 e/A 3 . Spin density of n-CiBN with n 
equals to (d) 8, (e) 10, and (f) 12. (g) and (h) are schematic 
energy diagram for C^BN. (g) Charge polarization caused by 
7r electron hybridization between C and BN. The blue arrow 
indicates the direction of electron transfer, (h) Competition 
between spin and charge polarization. Spin degenerated or- 
bit als are represented by long level lines, and spin orbit als are 
denoted by short level lines. Red and blue represent the spin 
up and spin down channel, respectively. 
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FIG. 3: (a) ~ (e) Configurations and density of states of a 
pair of B and N atoms substituted ZGNR. Fermi level is plot- 
ted with dotted line. Grey, blue, khaki and white balls denote 
carbon, nitrogen, boron, and hydrogen atoms, respectively. 



tory but it can serve as a proof of principles for achieving 
metal-free half metal. For the purpose of practical (mass) 
production of the ZGNR-based half metal, we propose an 
alternative material design as illustrated in Fig. [3] (left 
panels). Here, a pair of B and N atoms are doped in place 
of two C atoms in each supercell of ZGNR (dilute substi- 
tution). The location of the substitution is arbitrary. In 
Fig. [3] (right pannels), we show the calcualted DOS for 
five typical configurations of the BN substitution, three 
at the edge and two in the interior of the ZGNR. It can 
be seen that with the exception of the symmetric sub- 
stitution at the opposing edge [Fig. 09(c)], four out of 
five configurations exhibit spin polarization also. Some 
hybrid nanoribbons can be viewed as spin-selective semi- 
conductor [e.g. in the case of Fig. [31(b)] . 

The most interesting case is the configuration shown 
in Fig. EJe), where both B and N atoms are located in 
the interior of the ZGNR but not as a nearest neighbor. 
Appreciable DOS can be seen near the Fermi level in the 
spin-up channel whereas a finite band gap can be seen 
in the spin-down channel. The calculated band structure 
(left panel in Fig. [4j) shows that there is a very nar- 
row band gap in the spin-up channel (about a few tens 
of meV). Hence, this heterogeneous nanoribbon exibits 
half- metal or half semi- metal like behavior. This result 
suggests that half metallicity can be achieved by doping 
B and N atoms within the interior region of the ZGNR. 
The doping can be random as long as it is dilute to assure 
certain separation between the B and N atoms. To con- 
firm this design, in Fig. 0] (middle and right panels), we 
show the calculated band structures for two additional 
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FIG. 4: Band structures of randomly doped ZGNR with 
different width. Red for the spin-up channel, and black for the 
spin-down channel. Grey, blue, khaki and white balls denote 
carbon, nitrogen, boron, and hydrogen atoms, respectively. 

doping configruations in a wider ZGNR. Evidently, the 
band gap in one spin channel is narrower. This result 
suggests that electrons in this spin channel behave like 
metal (or semi- metal), especially when the the B and N 
atom is located near the two opposing edges, respectively 
(right panel in Fig. 2J. 

In conclusion, we have studied spin-dependent elec- 
tronic properties of heterogeneous C/BN nanoribbons us- 
ing density functional theory. We find that a ZGNR can 
be converted to half metal (or half semi-metal) either by 
inbedding an ordered BN row in ZGNR or by substitut- 
ing B and N atoms in the interior of ZGNR. In both 
cases, the tt orbital hybridization between C and B (or 
N) breaks the symmetry of the two edge states. Finally, 
we note that in both material designs no metal dopants 
are involved. This metal-free half-metallicity points to 
a new venue for making spintronic/electronic devices at 
low-cost. 
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